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Abstract 

Direct neutron capture on 62 Ni is calculated in the DWBA and the cross sections in the energy 
range relevant for s-process nucleosynthesis are given. It is confirmed that the thermal value of 
the capture cross section contains a subthreshold resonance contribution. Contrary to previous 
investigations it is found that the capture at higher energies is dominated by p- waves, thus leading 
to a considerably increased cross section at s-process energies and a modified energy dependence. 
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The energy range of about 10 to 50 keV is relevant for neutron captures in s-process 
nucleosynthesis in AGB stars and type II supernovae. A large number of efforts have been 
undertaken to provide the necessary information on nuclear cross sections and reaction rates 
for neutron captures on stable nuclides as well as long-lived isotopes in s-process branchings. 
Nevertheless, in a number of cases only thermal cross sections are known experimentally 
which have been extrapolated to the energy region of interest. This is often done assuming 
s-wave capture. However, there are two possible sources of error in such a treatment, even 
in absence of resonances in the relevant energy region. Firstly, the direct capture (DC) 
contribution may include higher partial waves which will become increasingly important 
with increasing energy. Secondly, the thermal cross section may be dominated by the tail of 
a subthreshold resonance. In the latter case, a simple s-wave extrapolation is not justified. 
For neutron capture on 62 Ni both effects are relevant. A calculation of the DC contribution 
can help to resolve both problems. 

The two most recent recommended values of the 30 keV Maxwellian Averaged Cross 
Section (MACS) are 35.5±4 mb and 12.5±4 mb, as quoted by Refs. fl and Q, respectively. 
A closer scrutiny shows that they are both based on the same experimental thermal neutron 
capture cross section but with different assumptions for the extrapolation to 30 keV 0. 
In the earlier extrapolation the thermal value was taken to contain DC only which was 
then extrapolated assuming s-wave capture. In the more recent re-evaluation it was realized 
that a narrow subthreshold resonance at -0.077 keV could have contributed to the thermal 
value of 14.5 barn. The new extrapolation was then performed by subtracting an estimated 
resonance contribution and assuming s-wave DC for the remaining cross section. 

This clearly shows the need for a more consistent treatment, including a more complete 
description of the DC process. In the following, the DC contribution to the thermal cross 
section is studied as well as the behavior in the range of several tens of keV. 

The DC calculation was performed in a way as adopted in many other investigations 

nh n 

(e.g. |4j,|5j), making use of the code TEDCA p, including El, Ml, and E2 transitions. See 
those references for more details on the DC formalism. The number of open parameters is 
considerably reduced by employing folding potentials for the optical potentials of target plus 
incident neutron and of the final neutron bound state. The folding potential is obtained by 
folding the projectile and target density distributions p a , pa with an energy- and density- 
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dependent effective nucleon-nucleon interaction v e g [?| 

V(R) = XVf(R) = A / / p a (ri)pA(r2)v c fi(E, p a , p Al s) df{drl , (1) 



with A being a potential strength parameter close to unity, and s = K + — Ti , where 
R is the separation of the centers of mass of the projectile and the target nucleus. For 
nucleons the density distribution is usually described by a delta function, simplifying the 
double integration to a single one. Often, the imaginary part of the potential is very small due 
to the small flux into other reaction channels and can be neglected in most cases involving 



neutron capture at low energies 



. In the present case this assumption is well justified 



because no other reaction channels are open in t 



re investigated energy range. The density 



D 



was derived from measured charge distributions 

For the potential describing the final bound state of the captured neutron the parameter 
Af is fixed by the known binding energy. The value of Ai in the entrance channel is left to be 
determined. This can be achieved by describing elastic neutron scattering data. In general, 
elastic scattering contains contributions from potential scattering and resonant scattering. 
Since a description of the non-resonant direct process is to be found, it is essential to use 
potential scattering cross sections or scattering lengths without resonant contributions. Here, 
A; is determined by reproducing the potential scattering length of 6.3±0.4 fm Q|. The 
experimental error can be included, leading to a range of permitted values for Ai bounded 
by the two extremes Ahi g h = 0.81552 and Ai ow = 0.81392. Thus, there remain no open 
parameters in this calculation which would have to be fitted to reaction data. 

Not only capture of neutrons into the ground state of 63 Ni is relevant but also transitions 
to excited states have to be included. Table 1 shows the considered final states and transi- 
tions, along with the adopted spectroscopic factors and the dominating partial waves. The 
level scheme of 63 Ni and the spectroscopic factors are taken from several 62 Ni(d,p) 63 Ni exper- 
iments, as usual. The spectroscopic factors are renormalization factors for each transition to 
a certain nuclear level, containing nuclear structure effects. The error on the spectroscopic 
factors is given as about 20%, therefore an additional uncertainty of 20% has been added to 
the error in A; arising from the scattering lengths. The relative importance of the contribu- 
tion of each transition depends on the reaction Q value, the possible partial waves (s, p, d, 
. . . -wave capture) and the spectroscopic factors. Therefore, nuclear states at high excitation 
energy do not contribute considerably, neither do states with small spectroscopic factors. 
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By far dominating are El s- and p-wave captures. Higher partial waves and multipolarities 
have been included in the calculation but as their contribution is lower by several ten orders 
of magnitude they are not explicitly considered in the discussion below. 

The total DC cross section is the sum of all transitions to the final states. Realizing that 
s-, p- and d- waves are dominant, a simple parameterization of the cross section can be given: 

X 



a 



YVE + ZE 3/2 , (2) 
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with the three terms referring to transitions with orbital angular momentum I =0, 1, 2. 
Including the errors derived above, the following range a\ ow < a < Ohigh is found: 

a high = 1.1 (^f^ + 3.85 x 10-^VE) , (3) 
a low = 0.9 (2^f^ + 3.525 x Kr 3 V£) , (4) 

with the center-of-mass energy E of the incident neutron given in keV and the resulting 
cross section a in barn. The factors before the brackets describe the uncertainty in the 
spectroscopic factors. In the above Eqs. OHU the d-wave terms are not shown because they 
are smaller by six orders of magnitude than the p-wave terms. With such a parameterization 
it is easy to disentangle the contributions from different partial waves. Contrary to previous 
estimates, p-wave capture dominates the cross section in the energy range relevant for s- 
process nucleosynthesis. Therefore, the cross section increases with energy after an initial 
decline, leading to a higher value at 30 keV. This contribution is entirely due to capture in 
the l/2 + state at E x = 2.955 MeV. This state is well established as it was seen in several 
experiments Q]. A number of 5/2 + states and a l/2 + state at higher excitation 

energy are also reached by p-wave capture but their contributions are small due to smaller 
Q values and smaller spectroscopic factors. 

Since potential scattering lengths were used, the derived potential should not include 
resonant contributions. In other words, the calculated thermal cross section of 5.21±0.76 
barn is the purely direct contribution to the measured cross section of 14.5±4 barn whereas 
the difference is due to the high-energy tail of the subthreshold resonance plus the low-energy 
tails of resonances at positive energies. The DC value also fits well into the thermal cross 
section systematics of the other Ni isotopes which do not show subthreshold resonances. 

It is quite simple to convert the above parameterization in Eq. (J2J) of the cross section 
a to one describing the MACS (a) = (av) /vt by multiplying the p-wave term by 1.5, the 
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6 Spin assignment from |12| . 
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FIG. 1: MACS in the energy range 1 < E c _ m4 < 100 keV. For the DC component, the upper and 
lower limit of the sum of all partial waves (total), and the high and low limit of the p-wave are 
shown. The results are compared to two values derived from thermal neutron capture data and 
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extrapolated to 30 keV, taken from []| and [2j|- A previous experimental value from [15|] is also 
shown. Finally, the contribution of resonant capture including several resonances and calculated 
from resonance parameters taken from ^(| is plotted (dot-dashed line). 

d-wave term by 3.75, and keeping the previous parameters, thus obtaining 

X 



DC 



E 



+ 1.5YVe + 3.75ZE 3/2 



(5) 



Thus, the DC MACS is considerably enhanced when including the p-wave capture. This 
leads to a MACS of 30.1 < (cr) DC < 40.2 mb at 30 keV. This is back to the original value 
given in Q| (but for different reasons), also consistent with older data 14, 2|. The resulting 
MACS up to 100 keV are shown in Fig. [TJ 

The resonant contribution can be estimated from the resonance parameters given in 
literature Q| and is shown as an additional line in Fig. 1. Although this may look close to 
a 1/v function at first glance, the behavior is not that simple because it is a sum of several 
resonances and the energy- dependence of the widths has been taken into account in the 
calculation performed with the multilevel i?-matrix code SAMMY For a determination 
of the total rate, it has to be added to the direct terms. As can be seen from Fig. 1, the 



resonant contribution is small (5.13 mb) at 30 keV. The fact that resonances are not very 
important at this energy is underlined by a comparison to statistical model calculations 
la ] which shows that the direct reaction mechanism is dominating. The level density in 
this nucleus is theoretically estimated to be close to the limit to allow application of the 
statistical model. Thus it is to be expected that any statistical model will even overestimate 
the actual statistical resonant cross section. Single, previously unresolved resonances could 
still contribute but due to the rather wide level spacing their effects might be limited because 
of the averaging taking place in the integration performed in the calculation of the MACS 
and the reaction rate. Nevertheless, a measurement in the relevant energy region is highly 
recommended. 

Neglecting d-wave capture and utilizing the above parameterizations, the resulting astro- 
physical reaction rate for the direct contribution in cm 3 mol -1 s _1 is given by 

N A (av) DC = 2M545 * 10? f X + 1.2926 x 10 2 Y T 9 ) , (6) 



with the reduced mass /i, Avogadro's number Na and the temperature Tg in 10 9 K. The 

constant term arises from s-wave capture, the temperature dependent term includes the 

p-waves. The thermal population of excited states in the target can be neglected at low 

temperatures. The rate (also the cross section and the MACS) is bounded by X[ ow = 

2.232 x 10~ 2 barn keV 1 / 2 , Fi ow = 3.1725 x 1(T 3 barn keV" 1 / 2 and A high = 2.985 x 1(T 2 

barn keV 1 / 2 , Ihigh — 4.235 x 10~ 3 barn keV -1 / 2 . This rate is a lower limit because resonant 

contributions are neglected. The rate obtained by fitting the resonant contribution shown 

in Fig. 1 has to be added to the above DC rate: 

. . / 0.0512461 1.63596 1/3 
N A M res = exp -3.80993 j^- + 21.4001T 9 1/3 - 3.29138T 9 

+0.3822399T 9 5/3 - 5.52698 In Tg) . (7) 



The fit function was chosen according to the standard format as specified in |18[. The 
resulting total rate Na (cv) = Na (<tv) dc + Na (w) vcs is valid in the temperature range 
0.001 < kT < 100.0 keV. 

The enhanced rate has some impact also on type II supernova nucleosynthesis. In a recent 



study of nucleosynthesis in type 



all stable nuclides up to Bi 
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I super novae, using an extended reaction network including 
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2l| . the smaller value of [2j was used. Although AGB 



stars are considered to be the main source of the s-process nuclides in the solar system (see 



e.g. 



22J), considerable s-processing also occurs in the early burning stages of massive stars 



producing the weak s-process component. (This component also provides seed nuclides for 
the 7-process which produces proton-rich nuclei by photo-disintegration.) In the supernova 
nucleosynthesis calculation an enhanced production of 62 Ni relative to the other Ni isotopes 
was found (see Fig. 2 in |20| and Figs. 2-6 in 21(). Using the previous value for 62 Ni(n,7) 63 Ni 
as given by [l| - being higher by about a factor of 3 - there is no such overproduction (see 
e -g- 3)- Using the rate proposed here in type II supernova nucleosynthesis calculations, a 
similar result can be expected because the new rate is - within errors - compatible with the 
previous value of 1[. 

Summarizing, the consistent treatment of the reaction 62 Ni(n,7) 63 Ni in the DC formal- 
ism confirmed the importance of a subthreshold contribution to the thermal cross section. 
However, it also showed that p-wave capture is the dominating contribution to the cross 
sections and rates in the energy range relevant for s-process nucleosynthesis, contrary to 
earlier assumptions. The new rate leads to 62 Ni abundances in type II supernovae consistent 
with the abundance levels of neighboring isotopes and curing a pronounced overproduction 
observed in some of the simulations. This underlines the fact that thermal neutron capture 
cross sections often cannot be simply extrapolated to the 30 keV region by assuming s-wave 
capture. On the other hand, the example of 62 Ni also shows that even thermal cross sections 
should be evaluated including a DC analyses to disentangle resonance effects and direct 
contributions. 

The remaining uncertainty in the cross sections can be reduced by resonance measure- 
ments in the relevant energy region to include possible effects of single resonances, by neutron 
scattering experiments with improved accuracy to reduce the error in the potential scatter- 
ing lengths, and by activation experiments within the relevant energy region which include 
direct and resonant contributions simultaneously but can currently only probe the MACS 
at a limited number of energies. In compilations such as [2| it would also be advantageous 
to identify those rates which are based only on the capture of thermal neutrons or are 
dominated by extrapolated terms. 
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